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Abstract 

We present a method, based on the positive weight next-to-leading-order matching formal- 
ism (POWHEG), to simulate photon production processes at next-to-leading-order (NLO). 
This technique is applied to the simulation of diphoton production in hadron-hadron col- 
lisions. The algorithm consistently combines the parton shower and NLO calculation, 
producing only positive weight events. The simulation includes both the photon fragmen- 
tation contribution and a full implementation of the truncated shower required to correctly 
describe soft emissions in an angular-ordered parton shower. 
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1 Introduction 

The production of photons via perturbative processes is very important for both the 
search for the Higgs boson and other new physics, via photon pair production, and 
for the study of QCD and experimental effects, in particular the jet-energy scale, 
in the production of a photon in association with a jet. To study these processes 
in detail in hadron-hadron collisions we need an accurate Monte Carlo simulation. 
In this paper we present a new approach for the simulation of these processes and 
illustrate it with the simulation of photon pair production. 

Monte Carlo event generators simulate events by combining fixed-order matrix el- 
ements, parton showers and hadronization models. The first programs used leading- 
order (LO) matrix elements, together with the parton shower approximation which 
describes soft and coUinear emission. Recently different approaches correcting the 
emission of high transverse momentum, px-, partons have been introducecQ. 

Different algorithms have been developed to provide a better description of the 
hardest emission, including the full next-to-leading-order cross section. In the ap- 
proach of Frixione and Webber (MC@NLO) [3ll[32], the parton shower approxima- 
tion is subtracted from the real emission contribution to the next-to-leading-order 
cross section and combined with the virtual correction. This method was success- 
fully applied to many different processes [33H10] . However, this approach has two 
drawbacks: it generates weights that are not positive definite and its implementation 
depends on the parton shower algorithm. 

These drawbacks have been addressed with a new method introduced by Na- 
son [1TP2] . the POsitive Weight Hardest Emission Generator (POWHEG) approach. 

^ See Ref. p] for a recent review of the older techniques [2Hl9] and techniques for improving the 
simulation of multiple hard QCD radiation [20ll30) . 
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1 Introduction 



This method generates positive weights and is implemented in a way that does not 
depend on the details of the parton shower algorithm. Nevertheless, the parton 
shower algorithm must have a well defined structure: a truncated shower simulating 
wide angle soft emission; followed by the emission with highest transverse momen- 
tum (pt)', followed by a vetoed shower simulating softer radiation. The hardest 
emission is generated by a Sudakov form factor that includes the full matrix ele- 
ment for real emission. The truncated shower generates emission at a higher scale (in 
the evolution variable of the parton shower), while the vetoed shower simulates ra- 
diation at a lower evolution scale than the one at which the hardest emission is 
generated. The POWHEG method has been successfully apphed to a wide range of 



These approaches have yet to be applied to processes involving the production 
of photons due to the complications which arise in the experimental measurement, 
simulation and calculation at higher orders in perturbation theory of these pro- 
cesses. Collider experiments do not measure inclusive photons because of the high 
background due to the production of photons in meson decays. Indeed, the in- 
clusive production rate of high px vr°, 77, u mesons is orders of magnitude bigger 
than for direct photon production. For this reason the experimental selection of 
direct photons requires the use of an isolation cut. Different criteria for the isolation 
of photons include: the cone approach [HSIET], the democratic approach [HB] and 
the smooth isolation procedure [69]. In fixed-order calculations this contribution 
is included using the measured photon fragmentation function, the probability of a 
parton fragmenting to produce a photon with a given fraction of the parent parton's 
momentum, whereas Monte Carlo simulations instead rely on the parton shower and 
hadronization models to simulate this contribution. This presents a problem in sim- 
ulating these processes at NLO where some of the singularities in the real emission 
processes are absorbed into photon fragmentation function in fixed-order calcula- 
tions. In this paper we will present a method for simulating these processes using 
the POWHEG approach which still relies on the parton shower and hadronization 
models to simulate the photon fragmentation contribution. This approach is simi- 
lar in its philosophy to the method of Ref. [TD] for combining leading-order matrix 
elements and the parton shower. 

We illustrate this approach with the simulation of diphoton (77) production. 
Diphoton production is important as it provides a large background for the discov- 
ery of the Higgs boson decaying into a pair of photons, for both the Tevatron [71] 
and LHC experiments [72l[73]. It is also an important background in new physics 
models, for example in heavy resonance models [71], models with extra spatial di- 
mensions [7S] and cascade decays of heavy new particles [7B]. Experimental mea- 
surements of 77 production have a long history in fixed-target [77H75] and collider 
experiments [801 - IM] . 

The theoretical understanding of diphoton production and precise measurements 
of the differential cross section are therefore not only important for the discovery of 
new phenomena but also as a check of the validity of the predictions of perturbative 

^ There has also been some work combining either many NLO matrix elements ^63j or the NLO 
matrix elements with subsequent emissions matched to leading-order matrix elements |64il65) with 
the parton shower. 
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quantum chromodynamics (pQCD) and soft-gluon resummation methods. 

The dominant production method for direct photon pairs is leading order qq 
scattering {qq — )■ 77), ahhough the formally next-to-next-to-leading-order, C(a^), 
gluon-gluon fusion {gg — )• 77) process via a quark-loop diagram [85] can be im- 
portant, and even comparable to the leading-order contribution at low diphoton 
mass (M^^) [8lj, due to the large gluon parton distribution function. 

The 0{as) corrections to the qq — )■ 77 process includes the qq — )■ 77(7, 
gq 77g and gq — )■ 775 processes and corresponding virtual corrections. More- 
over, the contribution where the final parton is coUinear to a photon is calculated in 
terms of the quark and gluon fragmentation function into a photon [851186] . Given 
the behaviour of the latter functions, ~ — , these terms contribute to the same 

' Ota ' 

order as qq — >■ 77. The QCD corrections to the process are well known in the lit- 
erature [5B| IS7H^ . Fixed-order Monte Carlo programs, such as JETPHOX and 
DIPHOX [23], provide simulation for direct photon production together with the 
implementation of isolation cuts. 

The present paper is organized as follows. In Sect. [2]we introduce the POWHEG 
formulae useful for the description of our approach and our treatment of the photon 
fragmentation contribution. The calculation of the leading-order kinematics with 
NLO accuracy in the POWHEG approach is discussed in Sect. [21 In Sect. H] we 
describe the procedure used to generate the hardest emission. We show our results 
in Sect. [5] and finally present our conclusions in Sect. [61 

2 The POWHEG method 

In the POWHEG approach the NLO differential cross section for a given N-body 
process is 

da = 5($B)d$B 
where is defined as 

B{<^b) = fi(<f b) + Vi<l>B) + j [R{^B, - C{^B, d$K, (2) 

B{(^b) is the leading-order contribution, $5 the N-body phase-space variables of the 
LO Born process whereas $ij are the radiative variables describing the phase space 
for the emission of an extra parton. The real contribution, R{^b, ^r), is the matrix 
element including the radiation of an additional parton multiplied by the relevant 
parton flux factors, and is regulated by subtracting the counter terms C{^b,^r) 
which contain the same singularities as R{^b,^r)- In practice the counter term 
is usually composed of a sum over a number of terms, $_r), each of which 

regulates one of the singularities in the matrix element using approaches of either 
Catani and Seymour (CS) ^94j or Frixione, Kunszt and Signer (FKS) [95], i.e. 
C($B,$ij) = The finite contribution V{^b) includes the virtual 

loop corrections and the counter terms integrated over the real emission variables, 
which cancel the singularities from the loop corrections, and the coUinear remnant 
from absorbing the initial-state singularities into the parton distribution functions. 



^«(0) + ^R))d^R 



5($ 



B 
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2 The POWHEG method 



The modified Sudakov form factor is defined in terms of tlie real emission matrix 
element 



^nipr) = exp 



(3) 



where kTi^B^^R) is equal to the transverse momentum of the emitted parton in 
the soft and collinear limits. 

The POWHEG method is based on two steps: the N-body configuration is gen- 
erated according to B{^b) and then the hardest emission is generated using the 
Sudakov form factor given in Eqn. [31 Since B{^b) is defined as the NLO differential 
cross section integrated over the radiative variables, the event weight will not be 
negative. 

If the parton shower algorithm is ordered in transverse momentum we would 
generate the hardest emission first and then evolve the + 1 parton final-state 
system using the shower forbidding any emission with transverse momentum higher 
than that of the hardest emission. On the contrary for shower simulations which are 
ordered in other variables, such as angular ordering in Herwig++ [T91I96] , the hardest 
emission is not necessarily the first one. For this reason the shower must be split into 
a truncated shower describing soft emission at higher evolution scales, the highest px 
emission and vetoed showers simulating emissions at lower evolution scales; however, 
constraints are imposed to guarantee that the transverse momentum of the emitted 
particles is smaller than the one corresponding to the hardest emission [^B2]. 

In order to use this procedure for processes involving photons where the real 
emission matrix elements contain both QCD singularities from the emission of soft 
and collinear gluons and QED singularities from the radiation of soft and collinear 
photons we need to make some modifications to the approach. We start by writing 
the real emission piece as 

^r) = Rqed{^B, ^r) + i?QCD($B, $i?), (4) 

where 

Rqeb{^b, ^r) = nP^T '^^^ ^^^^ 

l^j ^QED l^j -^QCD 

contains the collinear photon emission singularities and 

contains the singularities associated with QCD radiationjf] Here the counter terms 
have been split into those -Dqcd which regulate the singularities from QCD radiation 
and those -Dq^j-, which regulate the singularities due to photon radiation. 

We can regard the real QCD emission terms as part of the QCD corrections to 
the leading-order process, whereas the QED contributions are part of the photon 
fragmentation contribution coming from a leading-order process with one less photon 



In practice the counter terms can be negative in some regions and we choose to use their 
magnitude in this separation in order to ensure that the real contributions are positive. 
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and an extra parton. We therefore modify the next-to-leading-order cross section 
for processes with photon production giving 



da= ^B{^b) + V{^b) + 

+i?QED($B,$fl)d$Rd$B- 



d$R \ d$B 



(6) 



There should also be an additional non-perturbative contribution with the convo- 
lution of the photon fragmentation function and the leading-order process with one 
less photon and an extra parton. 

We can now write the cross section for photon production processes in the 
POWHEG approach in the same way as in Eqn.[T] 



d(T = i?($s)d$B 
+B\^'bWb 



Aqcd(O) + Aqcd(A:t($b, $«))d$R 



Aqed(O) + 



^qed(*'b,*'r) 



(7) 



B 



AQED(/i;r($B,$'i?))d$R 



where B{^b) is now defined as 



5(<I>b) = {B{'^b) + V{<^b) + 



d$i? > d$B 



and B'{^'b) is the leading-order contribution for the process with an extra par- 
ton and one less photon with and $'^2 being the corresponding Born and real 
emission phase-space variables. 

The Sudakov form factor for QCD radiation is 



Aqcd(pt) = exp 



- / d^n^^^^^^^^e{kri^B,^n)-PT) 



Bi^B) 

and the Sudakov form factor for QED radiation is 

i?QED($'i?,$'i?^ 



(9a) 



Aqed(pt) = exp 



d<l>' 



B'{(l>' 



B) 



9{kT{<l>'B,'^'R)-pT) 



(9b) 



Both the direct photon production and the non-perturbative fragmentation con- 
tribution are correctly included. The non-perturbative fragmentation contribution 
is simulated by the parton shower from the B'{^'b) contribution when there is no 
hard QED radiation. 

The POWHEG algorithm is implemented for photon production processes using 
the following procedure. 

• First select either a direct photon production or a fragmentation event using 
B{^b) and B'{^'b) and the competition method to correctly generate the 
relative contributions of the two different processes. 



For a direct photon production process: 
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3 Calculation of B{^b) 



— generate the hardest emission using the Sudakov form factor in Eqn.lQaj 

— directly hadronize non-radiative events; 

— map the radiative variables parameterizing the emission into the evolution 
scale, momentum fraction and azimuthal angle, {qh, Zh,(f)h), from which 
the parton shower would reconstruct identical momenta; 

— generate the iV-body configuration from B{^b) and evolve the radiating 
parton from the starting scale down to qh using the truncated shower; 

— insert a branching with parameters {qh, Zh, (ph) into the shower when the 
evolution scale reaches qh] 

— generate pt vetoed showers from all the external legs. 
• For a fragmentation contribution: 

— generate the hardest QED emission using the Sudakov form factor in 
Eqn.[9bl 

— directly shower and hadronize non-radiative events, forbidding any per- 
turbative QED radiation in the parton shower generating the 
non-perturbative fragmentation contribution; 

— for events with QED radiation map the radiative variables parameterizing 
the emission into the evolution scale, momentum fraction and azimuthal 
angle, {qh, Zh, (ph), from which the parton shower would reconstruct iden- 
tical momenta; 

— generate the A^-body configuration from B'{^'b) and evolve the radiating 
parton from the starting scale down to qh using the truncated shower, but 
allowing QCD radiation with pt greater than that of the hardest QED 
emission; 

— insert a branching with parameters {qh, Zh, (ph) into the shower when the 
evolution scale reaches qh, 

— generate the shower from all external legs forbidding QED radiation, but 
not QCD radiation, above the px of the hardest emission. 

This procedure now includes the QCD corrections to the leading-order direct pho- 
ton production process and both the perturbative QED corrections to the photon 
fragmentation contribution and the non-perturbative contribution are simulated by 
the parton shower. 

In the next two sections we will describe how we implement this approach in 
Herwig++ for photon pair production. 

3 Calculation of B{^b) 

In this section we describe the 0{as) corrections to diphoton production. At leading- 
order, 77-production is described by the Feynman diagram illustrated in Fig.[T] 
Next-to-leading-order contributions yield 0{as) corrections coming from qq — )■ 775', 
gq — )• 77g and gq — )■ 779, together with the corresponding virtual corrections, as 



3.1 Real emission contribution 
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+1 ^ 2 



Fig. 1: Diphoton production at leading-order. 

shown in Fig. El These subprocesses contain QED singularities, corresponding to 
configurations where the final-state parton becomes collinear to a photon, which 
do not cancel when summing up the real and the virtual pieces of the cross sec- 
tion. As described in the previous section they are formally absorbed into a quark 
/i^)) or gluon {G^g{z, fi'^)) fragmentation function into photons, which define 
the probability of finding a photon carrying longitudinal momentum fraction z in a 
quark or gluon jet at scale fi for a given factorization scheme. This QED singular 
component is called the Bremsstrahlung or single fragmentation contribution. In 
our approach it is treated separately and simulated by showering the gq — )■ 7g or 
gq 7g within the Monte Carlo algorithm, see Fig. [HI as described in the previous 
section. At next-to-leading-order the same configuration appears in any subprocess 
in which a quark (gluon) undergoes a cascade of successive collinear splittings ending 
up with a quark-photon (gluon-photon) splitting. These singularities are factorized 
to all orders in a^, according to the factorization theorem. When the fragmenta- 
tion scale IX is chosen higher than any other hadronic scale, i.e. ~ 1 GeV, these 
functions behave roughly as ^7-^^ and therefore they contribute at leading-order. 

For a full study at NLO accuracy, the 0{as) corrections to the Bremsstrahlung 
contribution need to be calculated. Moreover, these corrections in their turn yield 
the leading-order contribution of the double fragmentation type process; in the latter 
case, both photons result from the collinear fragmentation of a parton. However, 
these corrections are out of the scope of the present work and are not considered 
here. 

3.1 Real emission contribution 

In order to calculate the real emission contribution to B{^b) we need to specify both 
the radiative phase space, •I'/j, and the subtraction counter terms. We choose to use 
the dipole subtraction algorithm of Catani and Seymour ^94] to specify the counter 
terms and the associated definition of the real emission phase space as follows. 

In the centre-of-mass frame the incoming hadronic momenta are, P® and Pq, 
respectively for the hadrons traveling in the positive and negative z-directions. Sim- 
ilarly the momenta of the incoming partons in the Born process are = x^Pi^ and 
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3 Calculation of B{^b) 
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Fig. 2: Diphoton production at next-to- leading-order. In (a) the real and virtual 
Feynman diagrams contributing to the qq — )■ 77 subprocess are shown while 
in (b) the real diagrams for gq initiated process are given. 



pQ = XqPq, respectively. The momenta of the photons in the Born process are ki^2 
respectively. The corresponding momenta in the real emission process are and 
Pq for the incoming partons and ki^2,3 for the outgoing particles which are chosen 
such that ki^2 are the momenta of the photons and that of the radiated final-state 
parton. 

In the CS approach the real phase space depends on which parton is the emitter of 
the radiation and which the associated spectator defining the dipole [M]. When the 
parton with momentum p^ is the emitter and that with momenta pQ the spectator 
the full phase space is 



d$3 = d<l>Bd$/j = d$ 



{ki + k^f 



'B- 



167r2 2n 



X 



1 — X 



where the radiative phase space variables are 

(Pe +Pe) ■ h 



9{x{l—x)) 6{x—x^), 
(10) 



X = 1 — 



P(B-Pe 



P(B-Pe'' 



3.1 Real emission contribution 
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72 




Fig. 3: Bremsstrahlung contribution for diphoton production, 
is the azimuthal angle of the emitted particle around the ©-direction and 

xe[x(sA], V(se[0,l-x]. (12) 

In terms of these variables 



Pe = Pe^ 



(13a) 
(13b) 



It is useful to specify the momentum of the radiated parton in terms of its 
transverse momentum, pt, and rapidity, y, such that 



h = Pt (cosh y; cos 0©, sin 0©, sinh y) . 
Using the definition of x and f © we have 

h = V(sPe + (1 - X - V(s)P(B + q±, 



(14) 



(15) 



where q± is the component of the 4-momenta transverse to the beam direction. The 
on-shell condition, A;| = 0, gives 



- ?1 = Pt = 2pe ■ Pe(l - a; - ^^e)%. 
From Eqn. [15] and the definition of rapidity 



y = -In 
^ 2 



ftg -t- ftg 




I.E _ uz 
.'^3 '^3. 


2 



{1-x-v^)xq 

Vq)XXq 



(16) 



(17) 
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3 Calculation of B{^b) 



the CS variables are 



(18) 



ey 



This is sufficient to calculate the momentum of the radiated parton, however, 
rather than implementing the real emission variables in the Sudakov form factor in 
this way and then imposing the Oikxi^Bi ^r) — Pt) function it is easier to consider 
the real emission in terms of the transverse momentum, rapidity and azimuthal angle 
of the emitted parton. 

The Jacobian for this transformation is 



d{x, v) 



dprdy 



"iPT I 1 _ PT 



(l 4- Pt f^v\ ^ SX^Xq{1 — Vq 



The momenta of the photons in the real emission process can then be calculated 
from the Born momenta using 

k'^ = A''X r = l,2, (20) 

where the Lorentz transformation is 

with 

K = PcB + Pe - h = h + k2, (22a) 
K=P^+Pq. (22b) 

The condition = is compatible with the definition of x given in Eqn. [TTl The 
kinematic variables for the collinear direction are calculated in a similar way and 
they provide a radiative phase space as in Eqn. [TOl Moreover, given the x^ -H- Xq 
asymmetry of the rapidity in Eqn. [T71 it is [y]^ = — [y]^. In the rest of the paper 
we refer to the collinear direction as 6 = {0,©}, when both components need to 
be included. 

In addition to the real emission variables we need the dipole subtraction terms of 
Ref. [94j- In the following B{^b) and B'{^' b) are computed in terms of the reduced 
momenta defined in terms of the momenta for the real emission process in Ref. [91] . 
The QCD singularities from qq — )■ 77(7 are absorbed by the dipoles 

D'''' = = 87rC^a,(/i^)-^ - (1 + x)| 5($b), (23a) 



where the dipoles D'^^''^ denote the emitter i, emitted parton j and spectator k. 



3.1 Real emission contribution 
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The gq — > 775 subprocess involves the QCD dipoles 



D'"'" ^ D'^c^ = SnT^aM—- {1 - 2x{l - x)} B{^b] 



(24) 



In order to separate the QCD and QED emission we also need the QED dipoles 



-^97 — -^QED 



nil = n97^ 
— -^QED 



STrae 



87rae„ 



1 f 2 



-2 + z\b\^'b), 



B 



(25a) 
(25b) 



where 



kok, 



{k2 + /C3)p© ' 



{k2 + A;3)p© ' 



(26a) 
(26b) 



and Cq is the charge of the quark q in units of the electron charge. In this case, the 
radiative phase space is d$'^(^, z, (/)'). Similar dipoles are included for the gq — )> 77^ 
subprocess. We do not include perturbative QED radiation from the qq — > 751 
subprocess as it does not give a perturbative correction to G^g{z, ji'^). 

In practice we generate the real emission piece as a contribution from each of 
the incoming partons as 



/ 



(27) 











-^QED 




1 

-^qcdI 



d$ 



R- 



For the later generation of the Sudakov form factor it is useful to express the 
dipoles as 



QCD 



ZTT 



where / = {qg;qg; gq; gq}, 



Cqg — Cqg — Cp, 

^9Q ~ ~ "^Pi 



and 



(28) 



(29a) 
(29b) 



(30) 



where J = {q^yF, q^yl, q^yF, qyl}. 
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3 Calculation of B{^b) 



3.2 Virtual contribution and collinear remainders 

The finite piece of the virtual correction is 



dai 



27T 



-V{w)B{^b). 



where the finite contribution of 1(e) [9l] and the virtual correction [89] is 
V{w) = (3 + In^^ + In^l - w) + 31n(l - w)) + 

where Cq is the electric charge of quark q, and 

3 f 1 

F{w) = 21nw + 21n(l - w) + ^ ~ (Inw - ln(l - w)) 



w 



+ 2 



w 



1 — w 



In^u; + 2 



1 — w 



w 



with w = 1 + \, where s and t are the usual Mandelstam variables. 

s ' 

The collinear remainders are 

where the modified PDF i^ 

/-(xo,/iF) = £ ^ {/. M^) 



+ 



Xo 



fq -xfg{xo,fiF) 



Bix) 



+ fq (~'/^^) C'(^)} + fqiXO,fJ'F)D{Xo), 



fq and /„ are the quark and gluon PDFs respectively, and 



Up 



B{x) 
C{x) 



2x(l -x) + (x^ + (1 - x)2)ln^^^V^ 

fJjpX 

2 _j^g^(i-x)2 



1 — X 



2 , Q2(^_^)2 

1 — a; — ^ ma; — (1 + x)m ^ 



1 — X 



D{xo) 



-In + 21n(l - Xo)ln {\]+ 21n"(l - xq) + — - 5 



TT 



(31) 



(32) 



(33) 



(34) 



The combined contribution of the finite virtual term and collinear remnants is 



da 



y+coii 



2tt 



-V($i,)5($B) 



where 

with V{xo, A^f) 



V{<!>B) = V{w) + V{xo,fiF), 



(35) 

(36) 
(37) 
(38) 
(39) 

(40) 
(41) 



/{3;'o.mf) 



We write the modified PDF for the quark q, but a similar expression is vahd for an incoming 
antiquark q. 



3.3 Generation of the hard process 



13 



3.3 Generation of the hard process 

The next-to-leading-order simulation of photon pair production in Herwig++ uses 
the standard Herwig++ machinery to generate photon pair and photon plus jet 
production in competition. The B function is implemented as a reweighting of the 
leading-order matrix element as follows: 

1. the radiative variables and are transformed into a 
new set such that the radiative phase space is a unit volume; 

2. using the standard Herwig-(--|- leading-order matrix element generator, we gen- 
erate a leading-order configuration and provide the Born variables $ b with an 
associated weight B{^b)', 

3. the radiative variables are generated and B{^b) sampled in terms of the 
unit cube {x, v, 0), using the Auto-Compensating Divide-and-Conquer (ACDC) 
phase-space generator [57] : 

4. the leading-order configuration is accepted with a probability proportional to 
the integrand of Eqn.|8] evaluated at $_r}. 

4 The generation of the hardest emission 

Following the generation of the Born kinematics with next-to-leading-order accuracy 
the hardest QCD or QED emission must be generated according to Eqns.lQalor I9bt 
respectively depending on whether a direct or photon fragmentation contribution 
was selected. 

4.1 The hardest QED emission 

The hardest QED emission is generated by using the modified Sudakov form factor 
defined in Eqn. |9bl We generate it in terms of the variables $'^(xp, Zp, 0), with 



defined in [9l[T3], where 1], -2p G [0, 1] and the azimuthal angle G [0, 27r]. 

The invariant mass of the initial-final dipole = {pi — PkY — is preserved by 
the photon radiation. It is easiest to generate the hardest emission by introducing 
x± such that the transverse momentum of the emission relative to the direction of 



The Sudakov form factor can then be calculated in terms of (^'^{x±, Zp,(f)), such that 
the ^-function simply gives x± as integration limits and Eqn. [9b] becomes 



d$R = —dxpdzpd(f), 



(42) 



the partons in the Breit frame of the dipole is pr = ^x±, where 

2 4(1 Xp)(l Zp)Zp 



(43) 




(44) 
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4 The generation of the hardest emission 



where 



Z^j I^QEdI + l^j I^QCD 



the Jacobian, W, is 

>V = 4zp(l-^p)(l-Xp)^ (46) 

and yx™^^ is the maximum value for the transverse momentum. 

It is impossible to generate the hardest emission directly using Eqn.l44] instead 
we use an overestimate 

9ix±) = 4' (47) 
a; I 



of the integrand in Eqn.|44]so that 



- J ^d0d^,aj (48) 

can be easily integrated in {x±,x^^^}. This allows us to solve TZi = Aq^j^{x±) 
where TZi is a random number in [0, 1] to get the transverse momentum of a trial 
hard emission 

4(7^l) = 1 (49) 

This trial hard emission is then accepted or rejected using a probability given by the 
ratio of the true integrand to the overestimated value. If the emission is rejected the 
procedure is repeated with x^^^ set to the rejected x± value until the generated value 
is below the cut-off. This procedure, called the veto algorithm, correctly generates 
the hardest emission according to Eqn.HH [98j . 



4.2 The hardest QCD emission 

The hardest QCD emission is generated in terms of the variables ^r{pt, y, (p) defined 
in Sect. 13.11 Eqn. [Ha] then becomes 



^qcd(Pt) = exp - / dp'^d(l)dy 



W/^^ , (50) 



27r B 











-^QED 




-^qcdI 



where 

( 

'2% 

the Jacobian is 

= -r^. (52) 
1 - Vo 

where we mean to use f ® for / = {qg; gq; gq} and Vq for I = {qg}. 

As before we use the veto algorithm to generate the hardest QCD emission 
according to Eqn. [501 In this case we introduce the overestimate function 

9i{Pt) = — , (53) 
Pt 
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so that 

^qcd{Pt) = exp (- r ^d<pdyai) (54) 
V JpT Pt J 

is easily integrable in {pt^Pt^"^} ^^'^ — ^qcd(Pt) can be solved giving 

PT(7^l) = 7^^ (55) 

As before this trial hard emission is then accepted or rejected using a probability 
given by the ratio of the true integrand to the overestimated value. If the emission 
is rejected the procedure is repeated with -p^^^ set to the rejected value until the 
generated value is below the cut-off. 

5 Results 

Unlike the implementations of many other processes in the POWHEG formalism it 
is impossible to directly compare our results for any quantities directly with next- 
to-leading-order simulations in order to test the implementation due to the very 
different treatment of the photon fragmentation contribution. Instead we compare 
a simple observable, the rapidity of the photons, with the next-to-leading-order 
program DIPHOX [93] as a sanity check of our results not expecting exact agreement, 
although the PDFs and electroweak parameters were chosen to give exact agreement 
for the leading order qq — )■ 77 process. 

For proton-proton collisions at a centre-of-mass energy of 14 TeV, we used the 
following set of cuts on and rapidity of photons 

> 25 GeV, \y"'\ < 2.5, (56) 

together with a cut on the invariant mass of the 77-pair 

80 GeV < M^^ < 1500 GeV. (57) 

Moreover, we follow typical experimental selection cuts to isolate direct photons 
from the background: we require that the amount of total transverse energy, E)^^'^, 
released in the cone, centred around the photon direction in the rapidity and az- 
imuthal angle plane, is smaller than 15 GeV, i.e. 

{y - y^f + (0 - <P^f < R' (58) 
E^"'^ < 15 GeV, (59) 

where R = 0.4 is the radius of the cone. The PDFs are chosen to be the CTEQ6 
set [99]- The result is shown in Fig. HI The distributions from DIPHOX at NLO(red 
dashed line) and LO (red dash-dotted line), together with LO Herwig++ (dotted 
black line) and Herwig-I— |- with POWHEG corrections (solid black line) do not in- 
clude the gluon-gluon channel. At LO the Herwig++ and DIPHOX distributions are 
indistinguishable. At NLO they show a difference that is very small compared to 
the correction from LO to NLO, which means that the NLO curves are in reasonable 
agreement given the sizable contribution of the fragmentation contribution that is 
treated differently in the two approaches. 
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5 Results 
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Fig. 4: Rapidity of the 77-pair at NLO. The distribution from the Herwig++ parton 
shower with POWHEG correction (sohd black hne) is compared with NLO 
cross section from DIPHOX (dashed red hne). At LO the Herwig++ distri- 
bution is given by the dotted black line while the cross section from DIPHOX 
by the dash-dotted red line. 

In Fig.[5K we compare the results from Herwig++ with the data of Ref. [83] . 
a fixed next-to-leading-order calculation from DIPHOX (dotted magenta line) and 
RESBOS (dashed-dotted green hne) |100H104j . which performs an analytic resum- 
mation of the logarithmically enhanced contributions. Here and in the following the 
LO Herwig++ parton shower (red dashed line) includes the qq 77, qg — )■ 7jet 
and gg — ?■ 77 contribution. The implementation of POWHEG correction improves 
the description and this results in a distribution (solid blue line) that is in good 
agreement with the data. Here, as in the following, the NLO curve includes the 
gg — !■ 77 subprocess. In the lower frame, we plot the ratio MC/data and the yellow 
band gives the one sigma variation of data. All the plots comparing the results of 
Herwig++ with experimental results were made using the Rivet |105j package. 

It is of interest to study the transverse momentum of the 77-pair, because it 
is not infrared safe for p'Y — )■ 0. The qq — )■ 77 and gg — )■ 77 processes present a 
loss of balance between the corresponding real emission and virtual contribution, 
which results in large logarithms at every order in perturbation theory. In addition, 
the fragmentation components introduce an extra convolution that smears out this 
singularity. Since DIPHOX is based on a fixed, finite order calculation it is not 
suitable for the study of infrared sensitive observables and it fails in the description 
of these observables at low pJJ , as it is shown in Fig.|5b (dotted magenta line). 
Resummation for diphoton production in hadron-hadron collision has been provided 
at all orders in in Ref. |106j and implemented in RESBOS, as the corresponding 
distribution (dashed-dotted green line) shows in the same figure. The Herwig++ 
parton shower resums the effect of enhanced coUinear emission to all orders in 
in the leading-logarithmic (LL) approximation and results in a finite behaviour for 
pY — )• (red dashed line). However, the LO distribution does not correctly describe 
the data. In presence of POWHEG correction the distribution (solid blue line) stays 
finite at low pJJ and is in good agreement with the CDF data [83] . 

In addition, Herwig++ distributions, with and without POWHEG corrections. 




Hw++ POWHEG — 

Hw++ LO I 

--- DIPHOX NLO 
DIPHOX LO 
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Fig. 5: The (a) invariant mass and (b) transverse momentum of the 77-pair. The 
sohd blue hne shows the POWHEG approach, while the dashed red curve 
shows the result of the Herwig++ shower at LO. We show the NLO cross sec- 
tion provided by DIPHOX (magenta dotted line) and RESBOS (green dashed- 
dotted line). The data are from Ref. [83] and the curves are plotted with 
Rivet |105] . In the lower panel, the yellow band describes the one sigma 
variation of data. 
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(b) 80 GeV< M,,, <350 GeV 
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Fig. 6: Transverse momentum of the diphoton system for (a) 50 GeV < M^^ < 
80 GeV and (b) 80 GeV < M^^ < 350 GeV. The distribution for the 
POWHEG formalism (solid blue line) is plotted together with the distri- 
bution for the Herwig++ parton shower (dashed red line). The data are 
from Ref. [83] and the lower frame is as described in Fig.|5] 
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Fig. 7: Azimuthal angle between the photons for (a) 50 GeV < M^^ < 80 GeV and 
(b) 80 GeV < < 350 GeV. The sohd blue line shows the result for the 
Herwig++ shower with POWHEG corrections, while the red dashed line gives 
the result from the Herwig++ parton shower. The data are from Ref. [H3] 
and the lower frame is as described in Fig. [5] 
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(b) 80 GeV< M^^. <35o GeV 



-H- 



-H- 



-H- 



-H- 



2(7 1 I- 

"31^ hi I I I I I 



0.6 



0.7 

Icosfl'l 



Fig. 8: Polar scattering angle between the photons for two ranges of M^^: 50 GeV < 
M^^ < 80 GeV (a) and 80 GeV < M^^ < 350 GeV (b). The solid blue line 
describes the Herwig++ result with POWHEG corrections, the dashed red 
line does not include matrix element corrections. The data are from Ref. [83] 
and the lower frame is as described in Fig.O 
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are compared to the data of Ref. [HI]. In Fig.[6l we show the transverse momentum 
of the diphoton pair for two ranges of invariant mass of the 77-pair, M^^; in Fig.lH^ 
50 GeV < M^^ < 80 GeV and in Fig.Eb 80 GeV < M^^ < 350 GeV. For the 
same ranges of M^^ we plot the azimuthal angle distribution between the photons 
in Fig. [7^ and Fig.[7)D respectively and the polar angle between the photons in Fig. [8^ 
and Fig.[8]D. For all distributions we see that the LO Herwig++ ditributions (red 
dashed line) do not correctly describe the data. The POWHEG approach improves 
the simulation and provides a good description of DO data ^84j- 

6 Conclusion 

In the present work the POWHEG NLO matching scheme has been extended and 
applied to 77-production in hadron collisions. The QED singularities are not treated 
by including fragmentation functions but rather by simulating the LO cross section 
for the corresponding process and then showering it. The simulation contains a full 
treatment of the truncated shower which is needed to correctly generate radiation 
with transverse momentum that is smaller than the one of the hardest emission. 

The implementation of the process was tested by comparing the results with the 
fixed-order DIPHOX program which is in good agreement with the results of our 
approach for observables which are not sensitive to multiple QCD radiation. 

We find that without a correction to describe the hard QCD radiation there is a 
deficit of radiation in the simulation. The POWHEG approach overcomes this prob- 
lem and provides a good description of the data of Refs. . A remarkably good 
description is obtained for infrared sensitive observables, like the transverse momen- 
tum of the 77-pair at low pY , which demonstrates the resummation of logarithmic 
enhancement provided by the Herwig++ parton shower. 

This is the first NLO simulation of a process involving photons and provides an 
important new tool for the study of promt photon production. The simulation will 
be made available in a forthcoming version of the Herwig++ simulation package. 
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